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ABSTRACT 



We compute the synthetic 2-100 keV spectrum of 
an actively star forming galaxy. To this aim we use 
a luminosity function of point sources appropriate 
for starbursts (based on Chandra data), as well as 
the types of stellar sources and their correspond- 
ing spectra. Our estimates indicate that a Comp- 
ton spectral component - resulting from scattering of 
SN-accelerated electrons by ambient FIR and CMB 
photons - could possibly be detected, in deep IN- 
TEGRAL observations of nearby starburst galaxies, 
only if there is a break in the spectra of the brightest 
(L > 2 X 10'^® erg s~^) point sources. 
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1. INTRODUCTION 



Star formation (SF) leads to X-ray emission on var- 
ious spatial and temporal scales, including emission 
from O stars. X-ray binaries, supernovae and their 
remnants (SNRs), and galactic-scale emission from 
diffuse hot gas, and Compton scattering of FIR & 
CMB photons by SN-accelerated electrons. Inte- 
grated spectra of actively SF galaxies may show all 
these components as well as emission from an active 
nucleus (Rephaeli et al. 1991, 1995). 

Persic & Rephaeh (2002; hereafter, PR02) have 
quantitatively assessed the roles of the various X- 
ray emission mechanisms in SFGs. They have used 
an equilibrium stellar-population synthesis model of 
the Galactic population of X-ray binaries to deduce 
birthrates for interacting binaries; these, combined 
with estimates of the duration of the X-ray bright 
phase, have allowed PR02 to make realistic estimates 
of the relative (Galactic) abundances of high-mass 
and low-mass X-ray binaries (HMXBs and LMXBs, 
in which the optical companions are main-sequence 
stars with AI > 8 Mq and M < Mq, respectively). 
The abundance of SNRs (both Type II and la) 



was also consistently estimated. From the literature 
PR02 derived typical spectra for these source classes. 
The spectral properties and relative abundances of 
the various classes of stellar sources allowed PR02 to 
calculate the composite X-ray spectrum arising from 
a stellar population of Galactic composition. 

In addition to such "stellar" components, PR02 con- 
sidered two diffuse, non-stellar contributions to a 
starburst (SB) X-ray spectrum. 

(i) Thermal emission from the regions of interac- 
tion between SN-powered outgoing galactic wind and 
denser ambient ISM, mostly at energies <1 keV, and 
thus of little concern to us here. 

(a) Non-thermal emission due to Compton scatter- 
ing of the SN-accelerated relativistic radio-emitting 
electrons off the FIR and CMB radiation fields. The 
high SN rate in a SB is bound to yield high relativis- 
tic electron densities since SN shocks are known to 
be primary sites of cosmic ray acceleration. More- 
over, in a SB environment the mean FIR energy 
density can be ~ 10 times higher than that of the 
CMB. Compton scattering of a 100/i photon by a 
~1 CeV electron boosts the energy of the photon to 
~10 keV. A high density of electrons in a SB would 
obviously be reflected also in enhanced radio emis- 
sion. If the same electron population produces both 
the radio and hard X-ray emission, then both spectra 
are expected to have roughly the same power-law in- 
dex (Rephaeh 1979; Goldshmidt & Rephaeh (1995). 
Compton X-ray emission gauges the current SFR 
since the electron Compton-synchrotron lifetime is 
approximately <10^ yr (i.e. comparable to or shorter 
than a typical starburst duration). 

PR02 determined that the ^2-15 keV emission is 
dominated by LMXBs; this is not surprising, given 
that the synthetic spectrum was calibrated by that 
of a normal, quietly star forming galaxy like ours. 
Due to the long times between the formation and 
the onset of the X-ray phase (the optical companion 
is a subsolar main-sequence star) in LMXBs, their 
emission traces the higher, past Galactic SFR, and 
not the lower, current SFR. Thus in order to obtain 
an estimate of the current SFR, the spectral contri- 
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bution of LMXBs should not be included. As noted 
by PR02, at energies of approximately >30 keV (rel- 
evant for INTEGRAL) the SB spectrum is expected 
to be dominated by the Compton component. 

Recent observations of nearby galaxies - mostly with 
the Chandra satellite - have revealed several X-ray 
point sources (XPs) with 2-10 keV luminosities in the 
range 2 x 10^^ erg s~^ < L < 10^^ erg s~^ (hereafter 
dubbed Very Luminous Sources: VLXs), and a few 
Ultra-Luminous Sources (ULXs), with L > 10'^^ erg 
s"^ (Fabbiano & White 2003) ^ Although their na- 
ture is not yet known, VLXs and ULXs are observed 
to be directly linked to current star formation activ- 
ity (e.g., Zezas et al. 2002), and when identified, the 
optical counterpart of ULXs arc O stars (Liu ct al. 
2002), thus establishing a link between HMXBs and 
ULXs. The spectra of VLXs and ULXs are charac- 
terized both as power-law in the range of interest, 
with index F 2 and ^1.3, respectively. ^ 



2. CLASSES OF SOURCES OF X-RAY 
EMISSION 



Spectrum in 2-50 keV band: 

where e denotes energy (in keV), and with photon 
index 7 ~ 1.2, cutoff energy Cc ~ 20 keV, and e- 
folding energy ejr ~ 12 keV. 

d) Very Luminous X-Ray Sources (VLXs). 

Luminosity: 2 x 10^^ < L2-io/(ergs-i) < 10^^. 
Spectrum in 2-10 keV band: power-law with (photon) 
index F - 2. 

e) Ultra-Luminous X-Ray Sources (ULXs). 

Luminosity: 10^^ < L2-io/(ergs"^) < 5 x 10^^. 
Spectrum, in 2-10 keV band: power-law with (photon) 
index F ~ 1.3. 

f) Compton scattering. 

Luminosity: such that its flux at 10 keV is 0.1 of the 

total flux (sec discussion in Persic & Rephacli 2003). 
Spectrum: power-law with (photon) index /3 ~ 1.8. 



3. LUMINOSITY FUNCTION 



We briefly summarize the essential quantities needed 
to characterize each component of a SB spectrum: 



3 X 1032-33 erg s^^. 



a) O stars. 

Luminosity: ^2-10 

Spectrum: multi-T thermal {kT < 0.5 keV). 
(The contribution of O stars is very minor and will 
not be considered further here; see PR02.) 

b) Supernova Remnants (SNRs). 

Luminosity: 10^^ < L2-io/(ergs-i) < 10^^. 
Spectrum: thermal with kT ~ 2 keV, and Z ^ Zq. 

c) High-Mass X-Ray Binaries (HMXBs). 

Luminosity: < L2-io/(ergs-i) < 2 x lO^^. 

^These ranges correspond to Eddington luminosities for 
spherical accretion on, respectively, 2 and 8 Mq BHs. This is 
the range of BH masses that are thought to be attainable via 
ordinary stellar evolution. ULXs are called 'super-Eddington 
sources'. 

^Broadly speaking, the VLXs and ULXs could, in fact, 
describe different luminosity /spectral states of BH X-ray bi- 
naries (BHXBs) as a function of the mass infall rate (M) onto 
the BLl, while the other relevant parameters arc constrained 
within a much smaller range (e.g., BH mass M ~ 5 — IOMq, 
efficiency r) ~ 0.3 — 0.6). In fact, qualitatively speaking, 
the inverse correlation between luminosity and spectral in- 
dex (i.e., flatter spectra for higher luminosities) can be under- 
stood within the accretion model devised to explain the high 
X-ray luminosities of accretion-powered binary systems. The 
increase in luminosity is driven by an increase of M, which in 
turn means a piling up of material around the emission region 
that leads to a higher Compton optical depth, and hence to a 
lower value of V i.e. a flatter spectrum. Namely, if a source is 
embedded in a plasma cloud of optical depth t and temper- 
ature Te, the spectral index V of the high-frequency wing of 
the spectrum of the escaping radiation is 



3(r + 2/3)2feTc 



1/2 



The population of VLX and ULX sources in a galaxy 
is more abundant the more active is star formation 
in that galaxy. Speciflcally, if the point source in- 
tegral luminosity function (XPLF) is described as 
iV(> L) oc (which corresponds to a differen- 

tial XPLF of the type ^ oc i-(i+«)), based on a 
wealth of Chandra observations it has been found 
that XPLFs arc less steep in SFGs than in normal 
spirals and ellipticals. (For example, the two most 
nearby SB galaxies, have a ~ 0.5 (M82) and a. ~ 0.8 
(NGC253), while normal spirals have a ~ 1.2, and 
ellipticals have approximately a > 1.4; see Kilgard 
et al. 2002.) So the integrated point source luminos- 
ity of SFGs is dominated by the highest-L sources 
(Colbert et al. 2003). 



4. SYNTHETIC SPECTRUM 



Given an (observed) XPLF, with the "stellar" com- 
ponent clearly sorted by luminosity, we can formally 
compute the integrated " stellar" spectrum according 
to 



<^(e) = 



/ ^ i /(., L) dL 



dN(L) 



L dL 



(Sunyaev & Titarchuk 1980; Shapiro et al. 1976). 



/ dL 

where /(e, L) are the above spectral profiles, nor- 
malized to the energy band ei-e2 (here, 2-100 keV), 
where the luminosity L is defined. To this we add 
a diffuse Compton component. To be specific, we 
choose a power-law index a = 0.7 for the XPLF 
(i.e. roughly the average between the values observed 
for M82 and NGC253), and a Compton component 
which is 10% of the total fiux at 10 keV. It is impor- 
tant to realize that our extrapolation of the observed 
spectra to 100 keV is just for the sake of computa- 
tion. 
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Figure 1. The integrated 2-100 keV spectrum of a 
starburst. Shown in ascending order at 5 keV are 
the SNR, Compton, ULX, HMXB, and VLX com- 
ponents. The "stellar" components are normalized 
in energy flux in the 2-100 keV band. The Compton 
component has been normalized at 0. 1 of the total flux 
at 10 keV. The spectrum is photoelectrically absorbed 
through an HI column density N^i = 10^^ cm~^. 



then a ^^20 keV cutoff is expected also in the spec- 
tra of VLX and ULX sources. Hence, the knee of 
the synthetic spectrum will turn out more marked 
and the Compton component will possibly stand out 
more clearly (at energies of approximately >30 keV). 
INTEGRAL measurements can possibly determine 
whether there is a break in VLX/ULX spectra. 

On the other hand the Compton emission, in order to 
be discernible even in the absence of VLX/ULX spec- 
tral breaks, has to be stronger than assumed here. In 
our calculation we have assumed the Compton flux 
to be 0.1 of the total at 10 keV (which implies that 
it is -10% of the total 2-10 keV flux). This is a 
factor of ~2 higher than the upper limit set by Grif- 
fith et al. (2000) on the Compton flux fraction of 
M82 observed with Chandra (see Persic & Rephaeli 
2003). It is possible that in SB galaxies high SN 
rates and dust temperatures enhance the Compton 
component such that it dominates the spectrum at 
energies higher than ^-^30 keV. 



5. DISCUSSION 



The resulting 2-100 keV spectrum is shown in Fig.l. 
In the synthetic spectrum a conspicuous knee, due 
to the break in the HMXB spectrum, is immediately 
recognizable at 20 keV. If the ULX component, which 
turns out to lie very close to that of HMXBs for 
energies up to 20 keV, did turn out to have a similar 
break, then the knee in the synthetic spectrum would 
be sharper. 

An important related consideration is that the 
Compton component will be swamped by the VLX 
and ULX components and thus practically unde- 
tectable by INTEGRAL, unless both "stcUar" com- 
ponents exhibit a spectral break similar to that of 
HMXBs. It can be argued that the HMXB spectral 
break is related to the magnetic field strength, B, 
on the surface of the accreting neutron star: i) sta- 
tistically, the cutoff energy, Cc, is correlated with the 
observed cyclotron resonance energy, eqi by eg — 2 
as the latter depends on B so will the former, which 
will then provide a measure of B (a most helpful 
feature in those cases where a resonant scattering 
feature is not observed - White et al. 1995); and ii) 
physically, the spectral break could be due to radia- 
tive transfer effects in a strong magnetic field (Boldt 
et al. 1976). If the magnetic fields able to induce 
~20 keV spectral breaks in VLX-ULX spectra have 
strengths comparable to those observed in HMXBs ^, 

crude estimation suggests that for BH masses in the 
range 2 — 8 Mq , which can bo produced via ordinary evolution 
(i.e. through type II SN events) by progenitor stars with M ~ 
8 — 40 Mq , the conservation of magnetic flux implies that for 
_B ~ 10 — 100 G at the surface of the progenitor star, it will 
be S ~ 10^^ — 10^'' G at the Schwarzschild radius, nearly 
irrespective of the BH mass. (See Maeder & Meynet 1989 for 
the relevant quantities of the stellar progenitors) . These values 
are consistent with those deduced for the magnetized neutron 
stars which are members of HMXB systems (see Boldt et al. 
1976; White et al. 1995). 
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